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ABSTRACT 

The radio-quiet quasar BR1202-0725 {z = 4.695) is a remarkable source with a bright Northwest (NW) companion 
detected at submm and radio wavelengths but invisible in the optical. In the absence of amplification by gravitational 
lensing, BR1202-0725 would be the most luminous binary CO and far-IR source in the Universe. In this paper, we 
report observations with the IRAM Plateau de Bure interferometer of BR1202-0725 in the redshifted emission of the 
CO(5-4) and (7-6) lines, the [C I](^P2-^Pi) line, a high angular resolution (0.3" x 0.8") 1.3 mm map of the rest-frame, 
far-IR dust continuum, and a search for the CO(ll-lO) line. We compare these results with recent ALMA data in the 
[C II] line. Both the quasar host galaxy and its NW companion are spatially resolved in the molecular line emission 
and the dust continuum. The CO profile of the NW companion is very broad with a full width at half maximum of 
~ 1000 ± 130kms~^, compared to ~ 360 ± 40kms~^ for the quasar host galaxy to the Southeast (SE). The difference 
in linewidths and center velocities, and the absence of any lens candidate or arc-like structure in the field, at any 
wavelength, show that the obscured NW galaxy and the SE quasar host galaxy cannot be lensed images of the same 
object. Instead, we find morphological and kinematic evidence for sub-structures in both the NW and SE sources. We 
interpret these results as strong indications that the BR1202-0725 complex is a group of young, interacting, and highly 
active starburst galaxies. 

Key words, galaxies: high redshift - galaxies: individual: BR1202-0725 - techniques: interferometric - cosmology: 
observations 



1. Introduction 

' Submillimeter observations of galaxies and quasars at high 
[ redshift provide invaluable clues about the activity of 
. galaxy mergers leading to the formation of massive galaxies, 

■ and about the star-forming environment when the universe 
I was young. These clues come from the redshifted emission 

■ of dust heated by newly-formed stars and from molecular 

■ and atomic lines from the dense molecular gas in which 
. the stars are born. Since the early 1990s, massive galax- 

■ ies at high redshift have been observed with ground-based 
I millimeter and submm telescopes that can detect submm- 
. selected starburst galaxies and host galaxies of optically 

■ selected quasars. These detections are pushing more and 
more into the epoch of re-ionization, and are currently out 
to z = 7.1; (Venemans et al. 2012). 

One of the first high-z quasars that was detected and 
studied in detail in both thermal dust emission and in 
CO lines is the very luminous optically selected quasar 
BR1202-0725 at z = 4.7 (McMahon et al. 1994; Ohta et 
al. 1996; Omont et al. 1996; Benford et al. 1999; Carilh 
et al. 2002; lono et al. 2006; Riechers et al. 2006). The 
source is remarkable because it has two distinct compo- 
nents: a Southeast (SE) source associated with the opti- 
cally luminous quasar and an obscured Northwest (NW) 
source with no counterpart in the visible. The two submm 
dust continuum peaks, separated by 3.7" (linear distance 



24kpc)13 have nearly identical CO redshifts of z = 4.695 
(SE) and z = 4.693 (NW). The extreme infrared luminos- 
ity makes BR1202-0725 one of the brightest high-z sources 
{Liu — 3.7xlO^'*Lo, Leipski et al. 2010, confirming the 
original estimate by McMahon et al. 1994). The far-IR 
part of the luminosity is Lfir = 3.8xlO^^L0, with 1.2 and 
2.6x 10" Lq for the NW and SE galaxies (lono et al. 2006). 
The estimated star-formation rate in each component is 
> lOOOM0yr~^ and the combined molecular gas mass is 
- 10" M© (Omont et al. 1996; Riechers et al. 2006). 

The SE optical quasar host galaxy and the optically ob- 
scured NW source are both detected in the CO J = 1 — 
through 7—6 rotational lines (Ohta et al. 1996; Omont et al. 
1996; Carilh et al. 2002; Riechers et al. 2006) as weh as in 
the [C II] emission line, the main cooling line in the interstel- 
lar medium (lono et al. 2006; Wagg et al. 2012). At optical 
wavelengths, the SE component of BR1202-0725 appears as 
a quasar with at least two faint companion galaxies, seen 
in the optical (starlight) continuum (see our Fig. 6). The 
first of these optical continuum sources, found by Fontana 
et al. (1996) and Hu et al. (1996), is centered 2.6" north- 
west of the quasar, near to, but not coinciding with, the 
obscured NW submm source, and it appears to be part of a 

^ For linear sizes and luminosity distances we took the stan- 
dard cosmology with Ho — 71kms~^ Mpc~^, Q,m = 0.27, Qa = 
0.73 (Komatsu et al. 2011), and used the cosmological calculator 
by Wright (2006). 
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Fig. 1. Integrated CO (5-4) line maps, optimized to show total emission in the SE and NW galaxies separately. Source 
centroids and sizes are listed in Table 1. Both maps have natural weighting with a beam of 1.6" x 0.7" (PA 15°; lower-left 
box in left panel). Left: CO(5-4) integrated over a 444kms~^-wide band, centered on the SE galaxy's line center velocity. 
Contour steps are 0.24 mJy (Icr). The integrated CO(5-4) flux of the whole SE source, including the extension to the 
SW, is 1.6 Jykms"^. Right: CO(5-4) integrated over a 1067kms~"'^-wide band, centered on the NW galaxy's line center 
velocity. Contour steps are 0.15 mJy (Icr) The integrated CO(5-4) flux of the whole NW source, including the extension 
to the northwest, over the 1067 km s~^ band, is 2.6 Jykms"^ (see Table 2) 



filament of optical starlight continuum emission extending 
over 4" away from the quasar. Part of this optical contin- 
uum source is prominent in Lya emission at z =4.702 (Hu 
et al. 1996; Petitjean et al. 1996; Ohyama et al. 2004). The 
second, fainter, optical starlight continuum source is also 
a Lya emission galaxy, located 3" southwest of the quasar 
(Hu et al. 1997). 

Up to now, detailed study of the profiles of the molecular 
and carbon fine structure lines in BR1202-0725 have been 
hampered by lack of sufficient bandwidth to cover the very 
broad molecular emission lines of the NW galaxy (Omont 
et al. 1996; Carilli et al. 2002). To remedy this situation, we 
present in this paper new wide-bandwidth spectral-line ob- 
servations of CO and [Cl] emission in BR1202-0725, along 
with high-angular resolution observations of the 1.3 mm 
thermal dust continuum. These new data from the Plateau 
de Bure interferometer (PdBI) are compared with submm 
data obtained with ALMA that trace the [C II] emission line 
and the 0.9 mm thermal dust continuum emission (Wagg et 
al. 2012). Our observations improve significantly on previ- 
ous studies and shed new light on the kinematics of the 
SE and NW sources as well as on their structures. These 
data indicate that BR1202-0725 is a very complex group 
of merging starburst galaxies. 

2. Observations and data reduction 

The PdBI observed BR1202~0725 at 3 and 1.3 mm in 
February and March 2007 and at 2 mm in November 2011. 
In winter-spring 2007, the 6 antennas were in the extended 
A configuration with a longest baseline of 760 m. The re- 
ceivers were successively tuned to 101.190 and 222.478 GHz, 
close to the redshifted frequencies of the CO(5-4) and 
CO(ll-lO) lines. Due to the low declination of the source. 



the synthesized beam was elliptical with a FWHM of 
1.6" X 0.7" (PA 15°) at 3mm, and 0.85" x 0.26" (PA 15°) 
at 1.3 mm. 

The November 2011 observations were made in the in- 
termediate C configuration at 142.6 GHz, near the frequen- 
cies of CO(7-6) (redshifted to 141.637GHz) and [Cl](^P2- 
^Pi) (redshifted to 142.110GHz; hereafter, we refer to this 
neutral carbon line simply as C I). The beam was 3.0" x 1.8" 
(PA 18°) at 2 mm. 

The dual-polarization receivers operated in single side- 
band mode (LSB), with SSB receiver noise temperatures of 
40-50 K. The observations were done in dry weather, with 
system noise temperatures of ~140K at 2 and 3 mm and 
180 K at 1.3 mm, after correction for atmospheric absorp- 
tion and spillover losses. The IF bandwidth was 1 GHz in 
2007 (3000kms~i at 3mm, 1300 km s'^ at 1.3mm) and 
3.6 GHz in 2011 (7600 kms'^ at 2mm). In the raw spec- 
tra, channel spacings were 2.5 MHz in 2007 and 2.0 MHz in 
2011. 

Amplitude and phase were calibrated every 20 minutes 
on 3C 273, whose flux was compared on each observing day 
to the continuum of the ultra-compact HII region around 
the star MWC 349 (taken to be 1.04, 1.48, and 1.67 Jy at 
101, 141 and 222 GHz, respectively). 

Pointing and focus were checked every 40 minutes. The 
PdBI "seeing" parameter, which is not a limit on the radio 
resolution, but an indicator of the strength of the atmo- 
spheric phase fluctuations, was 0.33" to 0.49" at 3 and 2 mm 
and 0.23" to 0.26" at 1.3 mm. Short-term (< Imin) atmo- 
sphere phase fluctuations were corrected with the aid of 
3-channel monitoring receivers centered on the 22 GHz at- 
mospheric water line. Total on-source integration times, ex- 
cluding flagged data, were 13.2 h at 3 mm, 5h at 2 mm and 
6h at 1.3 mm. The data were calibrated with the GILDAS 
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Fig. 4. Position-velocity diagrams of CO(5-4) in BR1202-0725. Left: Dec. -velocity cut through the SE source. The main 
line has two peaks, at —70 and -t-110kms~^, separated by -fO.4" in dec. At negative velocities, a narrow-line CO tail 
extends > 3" southwest. This tail to the SW is near, but does not coincide with, the weak SW continuum source found 
by Wagg et al. (2012). Right: R. A. -velocity cuts through the NW source. Note the large linewidth, which includes the 
strong line core near — 200kms~^, and the weaker plateau, that may be another merger galaxy, at -1-400 km s""'^, and 
0.6" to the west or northwest. For both diagrams, velocity offsets are relative to 101.190 GHz (z = 4.6949), and contours 
start at 2 cr and increase in steps of 1-a (0.37mJy/beam). The beam is 1.6" x 0.7" (PA 15°); position offsets are relative 
to the map phase center at RA 12:05:23.11, Dec -07:42:32.10 (J2000). 



CLIC program; both polarizations were co-added and the 
maps were cleaned with the GILDAS MAPPING program. 

3. Results 

3.1. The 3-mm data: CO(5-4) line emission 

Figure 1 shows two different CO (5-4) integrated intensity 
maps, made separately to best show the SE and NW galax- 
ies, with a limited velocity range for the SE quasar host 
galaxy, and a much wider velocity range to show the total 
emission from the NW galaxy. These maps indicate that 
the CO in the NW galaxy is resolved to the north or north- 
west, in the direction of the beam minor axis (HPBW 0.7"). 
They also indicate a probable extension or tail somewhat 
south or slightly southwest of the SE galaxy. 

Figures 2, 3, and 5 show the CO(5-4) line data from 
the SE and NW galaxies in detail not previously seen. The 
most striking result is how dissimilar the spectra are. The 
CO profile of the SE galaxy, centered at redshift z — 4.6952, 
has a FWHM linewidth of AV = 360±40kms-i (Table 2), 
comparable to those of other high-z quasar host galaxies 
detected in CO. In the local universe, such a linewidth is 
typical of massive spirals seen edge-on, and in the mid- 
range of linewidths for merging galaxies like some of the 
low-z ultraluminous infrared galaxies (ULIRGs), and the 
high-z submillimeter galaxies (SMGs; see, e.g., Fig. 6 of 
BothweU et al. 2012). 

In contrast, the NW galaxy has an extraordinarily broad 
profile with AVfwhm — 1000 km s~^ and AVpwzp — 
1700 kms"^ (Fig. 4, right, and Table 2). These values ex- 
ceed the CO linewidths normally observed in quasars or 
galaxies, with some notable exceptions like the 1500 kms~^ 
-wide CO hne in the NGC 6240 merger (e.g., Feruglio et 
al. 2012), the high-z quasar-and-multiple-galaxy merger in 
SMM J02399-0136 (see Ivison et al. 2010 and references 
therein), and the z = 3.8 radio-galaxy merger in 4C41.17 
(De Breuck et al. 2005). The main CO line peak in the 



NW source is blue-shifted by — 200kms~^ relative to the 
CO centroid of the SE galaxy. This main CO peak domi- 
nates on the eastern side of the NW galaxy (Fig. 4, right). 
Toward the western side of the NW galaxy, there is a sec- 
ondary spectral feature that is red-shifted by -1-600 km s~^ 
relative to the main CO peak in the NW galaxy. These two 
spectral features at —200 and -1-400 km s""'^ suggest that 
the obscured NW submm object is in reality two merging 
galaxies, separated by about 0.6" (4.0 kpc). 

The redshift of the -1-400 km s^^ component happens 
to be the same as the z = 4.702 redshift of the Lya-l 
galaxy between the NW and SE galaxies (Hu et al. 1996; 
Fig. 6). Like the CO lines from the NW merger, the Lya 
line is also very broad with AVpvvzp — 1500 km s""'^ and 
AVfwhm ^ 1100 kms^^ (Ohyama et al. 2004). As in the 
CO in the NW and SE sources, there are also multiple veloc- 
ity components within the main Lya-1 galaxy, over a range 
of 500kms-i (Fontana et al. 1998; Ohyama et al. 2004). 
Hu et al. (1997) and Fontana et al. (1998) note that the 
C IV line, seen in the quasar, is absent in the Lya-1 galaxy, 
so the Lya is not reprocessed quasar radiation, but instead 
comes from star formation, at a rate of 15 to 54MQyr~^. 
The broad Lya linewidths may thus be due to superwind 
bubbles breaking out of the starburst and supernova re- 
gions in their parent galaxies. Because of these winds and 
because of optical depth effects, Lya redshifts can give very 
biased measurements of the systemic redshifts of galaxies, 
so the true redshift of the Lya-l object may differ from its 
measured Lya redshift by several hundred kms^^. 

As in the CO spectra (Fig. 2), the CO position- velocity 
diagram of the SE galaxy (Fig. 4, left) shows two main 
CO peaks separated by 0.4" to 0.5" in position, and by 
180 km s""'^ in velocity , that may indicate rotation of the 
molecular disk in the inner r = 2 kpc of the quasar host 
galaxy. There is also evidence of a weaker CO source or 
streamer extending 3" southwest of the SE galaxy. This 
tail is blueshifted by -100 km s""'^ relative to zero velocity 
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Fig. 2. CO(5-4) spectra across the SE galaxy, in steps 
of 0.4" in Dec. The tilted dashed line indicates how the 
two line components, at +110knis~^ in the north, and - 
TOkms"^ in the south, cause the line centroid to shift rela- 
tive to zero velocity (dashed vertical line through center of 
the spectra), in going from north to south (top to bottom). 
Position offsets in arcsec (Aa, A6, upper left of each panel), 
are relative to the map phase center at RA 12:05:23.11, Dec 
-07:42:32.10 (J2000). Velocity offsets (horizontal axes) run 
from -1200 to +1200 km s"! relative to 101.190GHz (z = 
4.6949). The beam is 1.6" x 0.7" at PA 15°. Channel widths 
are 59.3kms~^; la = 0.48 mJy; the data are not convolved 
with any smoothing function. 




R.A, 

Fig. 3. CO(5-4) maps in 178 kms~^-wide channels. The SE 
galaxy is strong at —44 and +133kms~^ (middle row, cen- 
ter and right panels), while the NW galaxy (upper right 
of each panel) has a main maximum at — 222kms~^ (left 
column, middle panel), and a weaker plateau extending to 
positive velocities, notably in the panel at 4-489 km s^^ 
(lower row, middle panel). Contour steps are 0.34 mJy; 
1(7 = 0.25 mJy. The cross in each box marks the position 
of the 1.3mm continuum from the NE galaxy (Table 1). 
Velocity offsets (upper left corner of each box) are relative 
to 101.190 GHz (z = 4.6949). The beam is 1.6" x 0.7" (black 
ellipse in bottom right panel). 



in Fig. 4, left, and it has a relatively narrow linewidth (100 
to 200kms~^). This weak CO tail to the southwest does 
not coincide, however, with the 0.9 mm continuum source 
reported by Wagg et al. (2012). Nor do we see, at the posi- 
tion of their weak SW continuum source, any CO emission 
at velocites corresponding to the edge of the bandwidth 
used by Wagg et al., where they had proposed a tentative 
detection of [C II] . 

The position- velocity diagram of the NW galaxy (Fig. 4, 
right) shows that the large linewidth in fact includes 
a strong core source centered at — 200kms^"'^, and a 
weaker plateau, that may be another merger galaxy, at 
-1-400 km s~^. This second line component, also seen in the 
channel maps of Fig. 3, is 0.6" west or northwest of the main 
source, in direction of the extension seen on the 1.3 mm dust 
continuum map. 

The 3-mm dust continuum is not visible in channel 
maps or in single-pixel spectra. From averaging the line- 
free channels on both ends of the spectra, we derive upper 
limits of 0.3 mJy (2a) for both the SE and NW sources. 
Spectral baseline fits to the spatially-integrated spectra 
yield a spatially-integrated, 3-mm continuum flux density 
of 0.2 ±0.1 mJy for the SE source, and are inconclusive for 
the NW source. 



3.2. The 2-mm data: CO(7-6) and CI lines 

Observations at 2 mm (142 GHz) show that the CO and [C I] 
lines are quite distinct in the SE source, where their spec- 
tral separation is larger than their linewidths. In the NW 
source however, due to their exceptionally broad widths, 
the CO(7~6) and CI lines are blended together (Fig. 5). 
For the NW source, because of this blending of the C0(7- 
6) and CI, and the two velocity components, we fixed 
the CO(7-6) linewidth to be the same as that of C0(5- 
4). The results for CO(7-6) and CI, from fitting a single 
Gaussian to each line are listed in Table 2. For the NW 
source, we show an alternative fit in Fig. 5, upper right and 
lower right, with two Gaussian velocity components in the 
NW source. In CO(5-4) (Fig. 5, upper right), this yields 
a main peak of 2.9±0.3mJy, at — 260±40kms~^, with a 
linewidth of 737 km s^^ (FWHM). There is a secondary 
peak at -K440±64kms"\ with a width of 268 km s"^ The 
two velocity components differ by 700±100kms~^, which 
agrees within the errors with the separation of the two ve- 
locity components in the channel maps of the NW source 
(Fig. 3), and in the position- velocity diagram (Fig. 4, right). 
The CO(7-6) line of the NW source (Fig. 5, lower right), is 
blended with the [C I] line, and a three- or four-component 
fit is needed. In the SE source, where the CO and [Cl] 
lines are well separated (Fig. 5, lower left), the [Cl] line 
has a width and center velocity consistent with those of 
the CO(7-6) and CO(5-4) lines (Table 2). 
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Fig. 5. CO(5-4), CO(7-6), and [Cl](3p2-3Pi) spectra of the SE and NW components of BR1202-0725. Histograms are 
spectra derived from fitting the interferometer visibihties by assuming elliptical Gaussian sources with sizes (FWHM) of 
1.0" X 0.6" (SE) and 0.8" x 0.6"(NW). Sohd curves show Gaussian spectral-line fits to the profiles, and dashed horizontal 
lines show continuum levels derived from these Gaussian spectral fits. The dotted histograms show alternative model-fit 
spectra for assumed point sources. The difference between model-point-source and model-extended-source spectra is more 
striking in CO(5-4) than in CO(7~6) because the larger beam in CO(7-6) more easily includes all the flux in each of the 
SE and NW sources. For the NW galaxy, two Gaussian velocity components are fit to the CO (5-4) and CO (7-6) lines, 
with a third Gaussian added to the CO(7-6) fit to include the CI line. Channel widths are 89kms^^ and 91kms^^ for 
CO(5-4) and (7-6) respectively, and velocity offsets (kms^^) are relative to CO(5-4) and (7-6) redshifted to z = 4.6952. 



For the 2 mm dust continuum, from both source maps 
and uv-fits, we obtain 2 mm continuum flux densities of 
1.3 ± 0.2 and 1.6 ± 0.2 mJy for the SE and NW sources. In 
the 3.0" X 1.8" beam at 2 mm, the dust continuum sources 
are unresolved. 



3.3. The 1.3mm data: dust continuum and search for 
CO(ll-lO) 

The high angular resolution of 0.26" x 0.85" (1.8x5.6kpc 
at z = 4.7) attained with the extended configuration of 
the PdBI at 1.3 mm enables us to study the thermal dust 
emission of BRl 202-0725 in detail. 

Figure 6 is a map of the 1.3 mm (222 GHz) continuum 
over the whole 1 GHz IF band. The band center was tuned 
to the redshifted frequency of CO(ll-lO). This line was 
not detected in the SE source, with a 3cr upper limit of 



5mJy/beam (Table 2). We see a 1.7cr hint of a line in the 
NW source, but more data are needed to confirm it. 

The dust continuum source positions and uv-fit sizes are 
listed in Table 1. The SE and NW sources are well separated 
in the 1.3 mm continuum. The peak positions of the dust 
emission are close to the peaks of CO (5-4) reported in this 
paper, and also to the CO (2-1) peaks (Carilli et al. 2002) 
and the non-thermal cm-radio continuum peaks (Momjian 
et al. 2005). 

The SE and NW sources are both resolved in the 1.3 mm 
continuum. Our model fits in the uv-plane yield equivalent 
circular Gaussian FWHM diameters of = 0.5" ± 0.1" for 
the SE dust source, and 0.4" ± 0.1" for the NW source. 
More elaborate fits are not justified at present; a good rule 
of thumb is that one needs S/N > 10 to measure, e.g., 
major and minor axes and position angles. The spatially- 
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Fig. 6. F5T image and 1.3 mm dust continuum map. Left: HST F814W image of the BR1202~0725 region (Hu et al. 
1996) showing the quasar and the two galaxies with Ly-a redshifts measured by Hu et al. (1996; 1997). In this filter, the 
two Ly-a galaxies appear only in the starlight continuum, not the Ly-a line. The field is 16.4" x 16.4". The smaller box 
shows the 5"x5" field of the 1.3 mm continuum map at the right. Right: Contours show the BR1202-0725 dust continuum 
at 1.3mm in steps of 0.26mJy/beam (l.Scr); negative contours are dashed, beam — 0.85" x 0.26". Next to the NW and 
SE galaxies separated by 3.7" are labels indicating the CO(5-4) velocity components discussed in the text, with offsets 
(kms~^) relative to 101.190GHz {z = 4.6949). The Ly-a-l object, near the NW galaxy, is indicated by a red ellipse. The 
cross in the red ellipse marks the optical continuum peak, 2.6" from the quasar (Hu et al. 1996; 1997). The Ly-a peak 
is near the southern part of the ellipse, 2.3" from the quasar (Hu et al. 1996; Fontana et al. 1998; Ohyama et al. 2004). 
In the visible and near-IR, there is also a string of starlight continuum emission extending over 4" in the same direction, 
but displaced by 0.6" from the Ly-a-l object (Hu et al. 1996; Fontana et al. 1998; Ohyama et al. 2004). In the SE galaxy 
(lower left) the central red dot indicates the X-ray source, coincident with the quasar (Evans et al. 2010). The blue dots 
show the VLA C0(2-l) peaks on either side of the quasar, and in the NW galaxy (Carilli et al. 2002). The light-blue 
diamonds show the extended non-thermal emission observed in VLBI at 1.4 GHz (Momjian et al. 2005). Source positions 
are listed in Table 1. 



integrated flux densities are 5.1±0.9 (SE) and 4.0±1.3 mJy 
(NW). 

The weaker dust source to the southwest of the quasar, 
reported by Wagg et al. (2012) at 340 GHz, is not seen in 
the 1.3 mm continuum map. Its 0.9 mm flux density of 1.9± 
0.3 mJy implies a 1.3 mm flux < 0.5 mJy, too low {2a) to 
be seen on our current map. 



4. Source properties 

In the Appendix to this paper, we review evidence that 
the SE and NW sources in BR1202-0725 are not strongly 
gravitationally lensed. All source properties given in the 
following sections therefore have no correction factors for 
lensing. 

4.1. The northwest source 

The above results show that the NW source is resolved in 
both the CO(5-4) line and the 1.3 mm continuum, as seen 
in Figs. 3, 5, and 6. The CO line proflles clearly change, 
between positions separated by 0.5", on either side of the 
1.3 mm continuum peak. The CO spectra and position- 
velocity diagrams across the NW source show that the 
sub-sources, or merging galaxies, have velocities that dif- 



fer by 600 km s ^, thereby explaining the unusually large 
linewidth of the overall CO proflle. 

4.2. The southeast source 

As noted above, the CO (5-4) spectra of the SE galaxy show 
two line peaks separated by 180kms~^, and 0.4" to 0.5" 
(Figs. 2 and 4). These two velocity peaks are probably the 
same as the two spatial peaks of the SE galaxy, already seen 
by Carilh et al. (2002) who used a VLA beam of 0.3" x 0.2" 
to resolve the C0(2-l) emission into two components, ori- 
ented SE-NW and separated by 0.3". These two compo- 
nents are on either side of the quasar (Fig. 6), and may 
be in opposite sides of a disk rotating around the nucleus. 
In the 1.3 mm dust emission, the quasar host galaxy also 
appears to have a weaker extension 0.5" to the northwest 
(projected distance 3.3 kpc; Fig. 6). 

From VLBI observations with the VLBA, the VLA, and 
the GBT, Momjian et al. (2005) claim that the SE 1.4 GHz 
nonthermal radio continuum peak is apparently 0.4" south 
of the quasar (Fig. 6, right, and Table 1), while the VLBI 
NW source coincides with the 1.3 mm NW dust peak. The 
SE displacement might be further evidence that the SE and 
NW sources have multiple components. The synchrotron 
emission sources seen in VLBI at 1.4 GHz are extended (di- 
ameters 0.3", or 2 kpc), with intrinsic brightness temper- 
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Table 1. Source positions and sizes. 



R.A. Dec. from Source 
Source 12''05'" -07°42' SE-1.3mm size Refs. 
(J2000) (J2000) continuum (FWHM) 



SE galaxy: 



SE 1.3 mm continuum 


23.129 


± 0.002 


32.74: 


± 0.03" 


n nil 
0.0 


0.52 ± 0.06 


1 


Chandra X-ray catalog 


23.130= 


± 0.02= 


32.74" 


± 0.3" 


0.0" 




2 


C0(2-l) SE (nw subpeak) 


23.122^ 


± 0.003= 


32.71" 


± 0.05" 


0.1" 


< 0.25" 


3 


C0(2-l) SE (se subpeak) 


23.134^ 


± 0.003= 


32.92" 


± 0.05" 


0.1" 


< 0.25" 


3 


CO(5-4) SE centroid 


23.130' 


± 0.002= 


32.94" 


± 0.06" 


0.2" 


0.5" ±0.1" 


1 


1.4 GHz nonthermal SE source 


23.1167 


= ± 0.04= 


33.109' 


" ± 0.3" 


0.42" 


0.30" X 0.17", 130° 


4 


NW galaxy: 
















NW 1.3 mm continuum 


22.978' 


± 0.002= 


29.79" 


± 0.03" 


3.7" 


0.26" ± 0.04" 


1 


CO(2-l) NW 


22.975= 


± 0.003= 


29.85" 


± 0.05" 


3.7" 


> 0.5" 


3 


CO(5-4) NW centroid 


22.971= 


± 0.004= 


29.67" 


± 0.06" 


3.8" 


0.7" ± 0.1" 


1 


1.4 GHz nonthermal NW source 


22.977= 


± 0.04= 


29.750' 


" ± 0.3" 


3.7" 


0.29" X 0.17", 115° 


4 



Position refs.: (1) This paper: our positions and sizes are from uv-fits to circular Gaussian models, and our 

absolute astromctric accuracy is 0.1"; 

(2) Evans et al. (2010); (3) Carilli et al. (2002); (4) Momjian et al. (2005). 



atures of 2 x 10^ K at their 8 GHz rest frequency. With a 
typical nonthermal spectral index of —2.7, they would 
have a rest-frame 1.4 GHz of 2 x 10^ K - about 500 times 
higher than in the central hundred parsecs of M82. This 
non-thermal radio brightness is comparable with those of 
compact starbursts in ULIRGs. The BR1202-0725 nonther- 
mal radio sources are therefore also likely to be the prod- 
ucts of extreme starbursts rather than active galactic nuclei 
(AGNs; Momjian et al. 2005). 

4.3. Molecular gas properties 

We estimated the molecular H2 gas densities and kinetic 
temperatures in the NW and SE sources by comparing 
the observed CO luminosities with single-component model 
brightness temperatures in the RADEX 1-D escape prob- 
ability program (Van der Tak et al. 2007). The data on 
CO (5-4), (7-6), and (11-10) are from this paper; the C0(2- 
1) intensities were derived from Carilli et al. (2002) (see 
Table 2). These are the only data sets that have both good 
sensitivity and clear separation of the SE and NW sources, 
and they thus provide the first opportunity to model the 
CO Spectral Energy Distribution (CO-SED) for the NW 
and SE sources separately (Fig. 7). 

We have no data from high density tracers like HCN and 
not yet enough spatial resolution for CO lines above C0(5- 
4). The kinetic temperatures and H2 densiticis arc therefore 
poorly defined. There are two main constraints: (1) in both 
sources, the CO(5-4) and (7 6) line luminosities, V , are 
nearly equal, and (2) in both sources, the C0(11 10) is 
very weak, and not detected, at least in the SE source, 
with a luminosity at least four to five times lower than the 
C0(5 4) luminosity (Table 2). 

If the gas temperature is in the range 40 to 50 K, i.e., 
typical of gas temperatures in ULIRGs, and close to the 
dust temperatures derived below, these constraints imply 
that the optical depths of the CO lines are high, at least 
up to the CO (7 6) line. Thanks to the collisional excita- 
tion and the resonant radiative trapping, the populations 



of levels J < 8 are then nearly thermalized (Te^ ~ T^), 
so that line brightness temperatures and luminosities L' 
remain constant. In higher rotational transitions, the line 
opacity decreases, and unless the H2 density is unrealisti- 
cally large (> lO^cm"^), this leads to a steep decrease in 
line brightness temperature. 

The key is the high optical depth in CO. The practical 
consequence for the BR1202-0725 sources is that in flux 
units, because of the weighting by v"^, the CO spectral-line 
energy distribution peaks around J=8 (Fig. 7). 

For comparison with escape probability models, we as- 
sumed the CO abundance was Xqq = 10"'', and we took 
the microturbulence to be lkms~^pc~^, both of which 
are standard values often assumed for Galactic molecular 
clouds. We took the background radiation field to be the 
cosmic background, at a temperature of 15.6 K at z=4.7, 
and we assumed gas kinetic temperatures of 43 K, the value 
estimated from the continuum dust spectra. In this zone of 
kinetic temperature, escape probability models can repro- 
duce the two constraints listed above for a nearly constant 
value of the product 



n(H2) 



= 1.2 X 10 



23 



(1) 



where n(H2) is the H2 density in cm~^, and Nqo/^v is the 
CO column density, in cm~^ per microturbulent linewidth 
At;, in kms~^. 

Realistic, optically-thick combinations would 
thus be n(H2) = 2xl03cm-3, and iVco/Au = 
6 X lOi^cm-2 (kms-i)-\ for both the SE and NW 
galaxies. The upper range would be n(H2) = 2xl0'*cm~^, 
and Nco/Av = 6 x 10^^ cm'^ (kms~^)~^ A solution in 
this range would be consistent with the gas mass (see next 
section), and with the dust becoming optically thick near 
rest-frame 100 fim. 

Lowering the CO abundance, or CO column density, or 
increasing the microturbulence would require higher densi- 
ties of molecular hydrogen to fit all the data. These simple 
models predict a C0(5 4) brightness temperature of 25 K, 
i.e., close to the Rayleigh- Jeans equivalent of (Te^ — T^g), 
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Table 2. Observed parameters of CO and atomic lines, and continuum fluxes. 



Line 


Reference 


Line 


Line 


Center 


Integrated 


Line luminosity 


Continuum 


and 


frequency 


peaf 


width^ 


velocity^ 


intensity*^ 




Scant ^ 


Source 


GHz 


mjy 


kms~^ 




Jy kms~^ 


10i°Kkms-*pc^ 


mJy 


CO (5-4) 


101.185 














oiii 




4.1 ± U.4 


o£;o _l_ 0'7 
ODO ± O ^ 


-1 1 -in 
— L( ± l( 


i.D ± U.Z 


o.o ± u.o 


U.z ± U.l 


JN W-main peak 








9fin -1- An 


Z.O ^ U.o 


7 Q _|_ 1 9 

/ .O ^ l.Z 


<, U.o yZj(j ) 


NW-extension 




0.9 ±0.1 


268 ± 50 


+440 ± 64 


0.3 ± 0.1 


0.9 ± 0.1 


< 0.3 {2a) 


NW-total 











2.6 ± 0.4 


8.7 ±0.1 


< 0.3 {2a) 


CO (7-6) 


141.637 














SE 




6.7 ± 0.5 


316 ± 40 


-27 ± 13 


2.3 ± 0.2 


3.9 ± 0.4 


1.3 ± 0.1 


NW-main peak 


— 


3.6 ±0.3 


737 ± 90 


-260 


2.8 ± 0.3 


4.8 ±0.5 


1.6 ±0.1 


NW-cxtonsion 


— 


1.2 ±0.1 


268 ± 50 


+440 


0.3 ±0.1 


0.9 ±0.1 


— 


NW-total 


— 


— 


— 


— 


3.1 ±0.4 


5.7 ±0.6 


1.6 ±0.1 


[C I](^P2--'Pi) 


142.109 














SE 


— 


1.7 ±0.4 


316 


— 


0.6 ±0.2 


0.9 ± 0.3 


1.3 ±0.1 


NW 


— 


1.8 ±0.3 


368 


— 


0.7 ±0.2 


1.2 ±0.4 


1.6 ±0.1 


CO(ll-lO) 


222.471 














SE 




< 3.8 (2(7) 






< 1.4 {2a) 


< 1.0 (2a) 


5.1 ±0.9 


NW 




2.7 ±1.7 






2.9 ±1.8 


2.0 ± 1.2 


4.0 ± 1.3 


[CII]'= 


334.96 














SE 




27 


328 ±6 


210 


9.6 ± 1.5 


3.0 ±0.5 


18 ±3 


NW 




20 


722 ± 12 


60 


14.7 + 2.2 


4.6 ±0.7 


19 ±3 


CO(2-l) 


40.479 














SE 




0.77 ±0.1 






0.23 ± 0.04 


4.8 ± 1.0 


< 0.16 {2a) 


NW* 




0.66 ±0.1 






0.39 ± 0.08 


8.1 ± 2.0 


< 0.16 {2a) 



Line peak fluxes, FWHMs, and velocities at the source peaks, from single-Gaussian spectral fits in velocity. 
Line intensities, L', and Scant are over 1.0" x 0.6" (SE) and 0.8" x 0.6" (NW); line luminosity L' from formula of Solomon et 
al. (1992). 

" [CII] data from Wagg et al. (2012). 

C0(2--l) flux in NW source (Carilli et al. 2002) corrected for limited VLA bandwidth by multiplying by 1.5, the ratio of the 
GO(5-4) flux in the full line width, to that in the limited VLA band of -525 to 175kms~^ 



the difference between the (Planck) excitation temperature 
and the cosmic background at the CO(5-4) rest frequency 
of 576 GHz. Figures 2 and 3 show that the peak C0(5 4) 
brightnesses are 2.0 and 1.5 mJy beam~^ in the SE and NW 
galaxies respectively. At 101.2 GHz, the 1.6" x 0.7" beam 
yields 107 K Jy"-'^, so in our beam, the observed peak bright- 
ness temperatures, corrected for (l+z), are 1.2 and 0.9 K 
for the SE and NW galaxies. If the source sizes at the ve- 
locities of the CO spectral peaks are about 0.3", or about 
the measured sizes of the dust sources, or about half the 
size of the entire line sources integrated over all velocities, 
then the true, deconvolved CO(5-4) peak line brightness 
temperatures, corrected for beam dilution, arc; 17 and 13 K 
for the SE and NW galaxies, or within a factor two of the 
predicted brightness temperature of 25 K from the simplest 
radiative transfer model. 



More complicated fits with multi-temperature and 
multi-density components would change these results, and 
meaningful values for Tk and n(H2) are hard to derive be- 
cause of the parameter degeneracy. In spite of the uncer- 
tainties, we can conclude, in agreement with previous stud- 
ies, that much of the molecular gas in both of the BR1202- 
0725 sources is dense and warm. 



4.4. Gas mass and dynamical mass 

We use the CO(l-O) line luminosity-to-(H2+He)-gas mass 
conversion factor derived for the center of ULIRGs by 
Downes & Solomon (1998), 

^gas/icOi_o = O.8M0(K km s'^ pc^)"! (2) 

to estimate the gas mass in the SE and NW galaxies. 

We took the C0(1 0) luminosities to be equal to the 
CO (5 4) luminosities. There; arc; two reasons to do this, 
observational and theoretical. Firstly, direct measiirements 
show that the CO(5-4) luminosities are equal to the C0(2- 
1) luminosities (Table 2), which means that in the CO 
levels J < 5, the CO is optically thick and thcrmalized. 
Furthermore, the total (SE+NW) CO(5-4) luminosity of 
1.3x10^^ Kkms^^ pc^ is not only the same as the total 
(SE+NW) CO (2-1) luminosity, but also the same within 
the errors as the total (SE+NW) CO(l-O) luminosity of 
1.0x10^^ K km s~^ pc^ obtained in the single-dish CO(l-O) 
observations by Riechers et al. (2006). 

Secondly, at z = 4.7, the CO(l-O) excitation tempera- 
ture and the gas kinetic temperature have to be warmer 
than the cosmic microwave background temperature of 
nearly 16 K for C0(1 0) to be in emission; otherwise C0(1- 
0) would be in absorption (up to now never observed in 



8 



p. Salome et al.: BR1202-0725: An Extreme Multiple Merger at z = 4.7 



E 



2 









1 1 1 1 

TOTAL 




' 1 










NW \ 






1 


1 


TV 
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Fig. 7. CO line peak flux densities vs. upper J level, for the 
SE (red lower curve) and NW (blue middle curve) galax- 
ies of BR1202-0725. The top black curve shows their sum. 
CO(l-O) data are from Riechers et al. (2006); C0(2-l) from 
Carilli et al. (2002); all other data are from this paper. 
The curves correspond to the radiative-transfer models de- 
scribed in the text. 



high-z sources), or have too low a contrast relative to the 
CMB to be detectable. In standard radiative transfer mod- 
els for gas warmer than 20 K (in order to have CO(l-O) in 
emission), and with niji-i) greater than 500 cm~^ (to col- 
lisionally excite the CO), and with high enough A^co (to 
make the CO opaque enough to be detectable), the bright- 
ness temperatures of CO(l-O) and 00(5-4) are the same. 
Hence at 2; = 4.7, at every place you detect CO(l-O), you 
will also detect CO (5-4) — with the same brightness. At 
z = 4.7, there is no extra "reservoir" of "cold" CO(l-O) 
that is currently detectable. 

We therefore take the CO(l-O) luminosities equal 
to the CO(5-4) luminosities, namely L' = (5.3 and 
7.8)xlO^°Kkms~^ pc^ for the SE and NW-main galaxies 
(Table 2), and thus derive Mgas = (4.2 and 6.2)x101°Mq 
for SE and NW-main, respectively. For both galaxies, the 
CO luminosities and consequently, the gas masses, are 
larger than in local ULIRGs, and comparable to the mean 
values of the SMGs listed by Bothwell et al. (2012; see their 
Fig. 6). For BR1202-07, these high CO luminosities and gas 
masses are probably partly the result of the sources' mul- 
tiplicity. 

Similarly, we can estimate the gas mass from the CI line, 
using the formulae of Weiss et al. (2003; 2005). Assuming 
the same excitation temperature as for CO, and a [CI/H2] 
ratio of 3xl0~^ by number (1.8xl0~^ by mass), we derive 
gas masses of (5.3 and 6.7) x 10^" Mq for the SE and NW- 
main galaxies, i.e., similar to the gas mass estimated from 
CO. 

Lower limits to the dynamical masses can be estimated 
from the CO linewidths and the source sizes. For the inner 
parts of the sources (the sizes derived from the fits to the 
1.3 mm continuum), we derive the dynamical masses of the 
SE and NW sources from 



Mdy„(< R) - —— ^232R- K'ot (3) 



where 1/G — 232, when Mdyn is in M©, R is in pc, and 
Viot is in kms~^ (e.g., Solomon et al. 1987). In principle, 
the rotation velocity could be estimated from the observed 
linewidths Al/pwRM, which would at face value imply true 
rotation velocities of 360 and 970 km for the SE and NW 
sources. In practice, these objects are dominated by the gas 
mass rather than the stellar mass, and have thick disks, 
with chaotic, merger-driven, and probably non-circular gas 
orbits. Above all, the turbulent velocities in the molecu- 
lar gas in advanced mergers and ULIRGs are very large 
, of the order of 100 to 150kms~^ (Downes & Solomon 
1998). Furthermore, the position- velocity diagram for the 
NW source (Fig. 4b) indicates that the NW CO source is 
multiple; the strong core of the CO line, with a width of 
~ 400kms~^ corresponds to the main NW galaxy and its 
1.3 mm dust continuum, while the plateau in the CO emis- 
sion to higher velocities may be the other galaxy in the 
NW merger, possibly the weaker satellite extension seen 
on the dust continuum map. If we take the 400kms~^ CO 
linewidth associated with the 0.6"-diameter CO source, and 
correct downward for turbulent broadening, the true rota- 
tion velocity would be 300 to 350kms~^, comparable to a 
rotation velocity we might deduce for the SE galaxy. Hence 
the dynamical mass within a radius of 2kpc (^ 0.3") in 
both galaxies would be of the order of 6 x 10^° M© for each 
of the SE and NW sources. 

Comparison of these dynamical masses of 6 x 10^° Mq 
with the gas masses of 4.2 x 10^° and 6.2 x 10^° Mq indicates 
that in the inner parts of the SE and NW galaxies, most 
of the mass is in the form of molecular gas. It also means 
that it would be difhcult in these galaxies to have a CO- 
to-H2 mass conversion factor much larger than the value 
for ULIRGs; otherwise, the gas mass within the central r 
— 2 kpc would greatly exceed the dynamical mass. 

4.5. Dust mass and infrared luminosity 

The dust properties in BR1202-0725 have been discussed 
by Bcnford et al. (1999), lono et al. (2006), Leipski et al. 
(2010) and Wagg et al. (2012). To these studies, the resuhs 
in this paper now add new measurements at 1.35, 2, and 
3 mm that update previous data from Omont et al. (1996) 
and Guilloteau (2001). For the NW and SE sources, the 
updated dust continuum fluxes are given in Table 2. 

For an estimate of the dust mass, we assumed a source 
size of 0.5", for consistency with the uv-fit sizes derived 
from the observed 1.3 mm continuum sources, and a dust 
mass opacity coefhcient of K{i'r) = 10 • i'^, where i^r is the 
rest-frame (emitted) frequency in THz, k is in cm^^ gm~^ 
of dust, and we took the index /3 = 2. This absorption co- 
efficient agrees with extensive literature estimates for very 
dense molecular clouds with coagulated dust particles (e.g., 
Kriigel & Siebenmorgen 1994; Ossenkopf & Henning 1994). 
Because our observing frequencies correspond to rest fre- 
quencies well into the far-IR, where the large quantity of 
dust starts to become opaque, it is important in these es- 
timates of dust mass not to use the standard formula for 
optically-thin dust. Instead, one must do a proper fit of the 
dust source function minus the background source function 
(see eq.(l) of Weiss et al. 2007). Because at z = 4.7, the 
cosmic background temperature is 16 K, it is important 
to subtract the background from the dust source function, 
as it reduces the source-to-background contrast for the ob- 
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Fig. 8. Dust continuum spectra for NE and SW galaxies 
in BR1202-0725. Flux densities from 100 to 340 GHz are 
from the mm/submm interferometer data, as these are the 
only observations that spatially separate the NW and SE 
galaxies in the cool dust continuum: 101, 142, and 222 GHz 
from this paper, 340 GHz from Wagg et al. (2012). Because 
the mm/submm dust continuum fluxes of the two galaxies 
are nearly the same, we plot the NW galaxy fluxes (blue 
lower curve) divided by five, for better readbility. Data 
points at 350, 160, 70, and 24 /im are from Benford et al. 
(1999), Leipski et al. (2010), and Hines et al. (2006). The 
red curves show model dust components for the SE galaxy: 
a starburst-heated 50 K component and two quasar-heated 
warm (130 K) and hot (410 K) dust components. The up- 
permost dotted curve is the total emission. 



served dust flux. Applying this method to the mm/submm 
interferometric data only, which are the only observations 
that spatially separate the dust continuum of the NW and 
SE galaxies, we can fit the observed data (Fig. 8) with 
a dust temperatures of 40 to 50 K in both sources. The 
luminosities of these "cool" , "starburst" dust components 
alone, from integrating over the range 20 to 1000 /xm (rest- 
frame) under the modified Planck functions in our models, 
are 6.3 x IO^^l© (SE) and 1.2 x 10" Lq (NW), compara- 
ble with values derived in previous studies. These values 
yield dust masses of 4 x 10® M© in both the NW and SE 
sources. 

The flux densities from 350 to 24/j,m (rest-frame 60 to 
4.2 /xm), follow the usual quasar-template monotonic de- 
crease (a flat curve when plotted as i^5y), due to quasar- 
heated dust from a continuous distribution of dust layers 
of increasing density and decreasing mass and size as one 
moves inward toward the quasar. Following the model of 
Weiss et al. (2003) for the Cloverleaf quasar, we approx- 



imate this continuous distribution with just two compo- 
nents, a "warm" (130 K) layer of size 0.11", and a "hot" 
(410 K) layer of size 0.014". This approximation (Fig. 8) 
shows that nearly all the rest-frame mid-IR and near-IR 
flux comes only from the SE galaxy (the quasar); the mid- 
to-near-IR contribution of the NW galaxy is negligible. The 
350 /im flux measured by Benford et al. (1999; rest frame 
60 /im), comprises not only the fluxes of the SE and NW 
"starburst" components, but also a substantial contribution 
(30 to 50 percent) from the quasar warm dust component. 
The Herschel FAGS fluxes at 160 and 70 /im (Leipski et al. 
2010; rest frame 28 and 12 /im) arc entirely from the quasar- 
heated "warm" and "hot" dust components; the contribu- 
tions from the SE and NW "starburst" components are 
negligible. The Spitzer MIPS flux at 24 /tm (Hines et al. 
2006; rest-frame 4.2 /tm) is mainly from the quasar-heated 
"hot" dust component. There is no contribution from the 
SE and NW "starburst" components. 

4.6. Star-forming activity and merging 

From the far-infrared luminosities derived above for the 
"cool", "starburst" dust components, one can estimate the 
star formation rates following Kcnnicutt (1998): 

SFR [Mo/yr] = 1.7 X 10-i°Lfir [Lq] (4) 

which yields 1000 and 2000 Mgyr"! for the SE and NW 
components of BRl 202-0725. The Kennicutt SFR relation 
is for a Salpeter Initial Mass Function (IMF). Adopting 
other IMFs, such as done by Kroupa (2001) or Ghabrier 
(2005), would reduce these SFRs by about a factor of two. 
Even with a more top-heavy IMF however, it is clear that 
both the SE and NW galaxies have very intense starbursts, 
in multiple active regions. Depending on the IMF, the min- 
imum times to consume all the molecular gas are 40 to 
80Myr (SE) and 30 to 60Myr (NW), but longer if there 
are several, spaced-out, million-year bursts. These times are 
comparable to the ~40 Myr-consumption times in typical 
submm galaxies (SMGs) at z - 2 (Tacconi et al. 2008). The 
derived gas masses in BRl 202-0725 are also about the same 
as in typical SMGs at z '-^ 2. It is possible, that the quasar 
in the SE galaxy, and possibly a radio-quiet, dust-obscured 
AGN in the NW galaxy, contribute to the dust heating, 
and the star formation rate is overestimated. The other 
possibility is that we are observing an advanced multiple 
merger in BR1202-0725 close to the peak in its starburst 
luminosity curve, as shown by the wealth of components 
in the BR1202-0725 merger, including possibly two sources 
within each of the SE and NW galaxies, the active quasar, 
and the companion Lyman a protogalaxies. This wealth of 
activity may at least in part be a signature of the evolution 
of massive, high-redshift galaxies. Similar complexity will 
probably be discovered in the SMGs a,t z ^ 2, once they 
have been studied in comparable detail to BR1202-0725. 

Rather than a merger of just two gas-rich galaxies, fa- 
miliar from many ULIRGs at low redshift, the complexity 
of the BR1202-0725 interaction, with multiple galaxies or 
pre-galactic building blocks, may itself be a characteris- 
tic of the assembly of massive galaxies at high redshifts. 
Another high-redshift example of a complex merger is the 
z = 3.8 merger in the radio galaxy 4G41.17, where there 
are at least two CO sources, with a total GO(4-3) line 
width of 1000 kms~^ around the powerful AGN source of 
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the extended radio jets (De Breuck et al. 2005), as well 
as numerous near-IR companions (Graham et al. 1994), 
that appear like raisins in a pudding contracting to form a 
very massive galaxy. At least one of the companion regions, 
4C41.17-South, has been shown to have the same redshift 
as the main AGN galaxy, 4C41.17-North (Van Breugel et 
al. 1999). Another relevant example is the z=2.8 multiple- 
galaxy merger in SMM J02399-0136 (see Ivison et al. 2010 
and references therein), which also shows 1000 km s""'^- wide 
CO lines coming from the merger of a far-IR luminous but 
highly obscured starburst (their L2SW galaxy, coincident 
with a Lya cloud), a BAL quasar host galaxy (their LI 
galaxy) and two other components in the field (the L2 and 
LINW galaxies). All four of these objects arc at the same 
redshift (Ivison et al. 2010), and must be kinematically as- 
sociated, not just seen in projection. As in BR1202-0725, 
the ensemble of merger components in SMM J02399-0136 
form a vast and complex nursery, not only for new stars, 
but also for the formation of the future massive galaxy it- 
self. 

In the new data, the NW galaxy and the SE quasar 
host galaxy each break up into sub-components. Based on 
the present work, the BRl 202-0725 group consists of at 
least four distinct sources (six when counting the faint Lya 
companions) . The proximity and similar redshifts of these 
galaxies or pre-galactic objects suggest that the SE quasar 
host galaxy, the obscured NW submm galaxy, and the Lya 
objects are all in a gravitationally-bound group of galaxies 
that have decoupled from the cosmological expansion, and 
are probably merging together. All of these objects thus 
have exactly the same cosmological redshift, and the slight 
differences in their measured redshifts are due to true kine- 
matic motions due to their mutual interaction within their 
"local group." 

In the standard cosmology, a z = 4.7 galaxy is at a co- 
moving (tape-measure) distance from us of 7.6 Gpc, and 
the lookback time is 12.4 Gy. When the galaxy emitted 
the light we now receive, its recession velocity was 2.3 c, 
the Hubble parameter, H{z = 4.7) was 505kms~^ Mpc~^, 
the only relevant component of the universe was matter, 
with fljii — 0.98, and dark energy was irrelevant, with Q,\ 
= 0.014 (see e.g., formulae in Ryden 2002; Liddle 2003). 
In the usual interpretation, the BR1202-0725 group is in 
the potential well of a overdensity of matter relative to the 
mean density of the universe, that, long before the z = 4.7 
epoch, had decoupled from the cosmic expansion, so the 
galaxies in the group have only true kinematic relative mo- 
tions governed by the total mass in the overdensity region. 
The greater neighborhood of the BR1202-0725 overdensity 
region as a whole may be similar to the numerous "redshift 
sheets" seen at high-2; in the Hubble Deep Field (Cohen 
et al. 2000), or to the recently- found galaxy groups at z 
= 5.2 (Capak et al. 2011; Walter et al. 2012). In interpre- 
tations of the same-redshift groups, (e.g., Cohen, Hogg, & 
Blandford 2000, their sect. 6), these regions extend over 
several Mpc and have virial masses of a few times 10^"^ to 
10^* Mq, resulting in true kinematic velocity dispersions 
(not cosmological redshifts) of a few hundred kms"-'^. Our 
measured line-of-sight velocity differences of the different 
BR1202-0725 objects are of this order, and due to the close 
approaches (projected separations 3 to 25 kpc) in the indi- 
vidual mergers, the relative velocities can be even higher. 

In other respects, the BR1202-0725 objects agree 
with the general picture of high-redshift protogalaxies. 



Numerous deep-field morphological studies show that at 
high- 2;, there are no spirals (e.g.. Driver et al. 1998; 
Abraham & van den Bergh 2002). The galaxies are gas- 
rich, puffcd-up disks or spheroids, and much smaller than 
modern-day galaxies. The data of, e.g., Oesch et al. (2010) 
show that the half-light radii decrease as (1+^) to the 1.1 
to 1.3 power. At z = 4.7, typical galaxy radii are 700 
to 1400 pc, depending on the mass. The high end of this 
range is close to sizes measured in the present paper for the 
BR1202-0725 SE and NW galaxies. 



5. Conclusions 

Our high-resolution millimeter observations of the BRl 202- 

0725 galaxy complex at z = 4.7, in both line and contin- 
uum emission, reveal new aspects of this extreme multi- 
ple merger. Morphological and dynamical evidence together 
with the lack of any lens candidate or arc-like structure in 
the field, at any wavelength, definitely rule out the hypoth- 
esis that the obscured NW galaxy and the SE quasar host 
galaxy are lenscd multiple images of the same background 
object. Instead, the present data indicate that the images 
of BR1202-0725 directly show a group of merging and very 
active starburst galaxies (Fig. 6, right). 

The SE galaxy merger contains at least two sources. The 
main SE core source coincides with the optical quasar's po- 
sition and shows internal kinematics compatible with rota- 
tion and a CO line profile (FWZP= 700 km/s) typical of a 
massive disk galaxy. To the southwest, the CO(5-4) line is 
narrower, with a linewidth of 200kms~^, and blueshifted 
by — 180kms~^ relative to the main SE velocity peak. This 
fainter CO does not coincide in position with the weaker 
SW continuum source reported in the ALMA data (Wagg 
et al. 2012), which may be a third, fainter, galaxy in the 
merger. 

The NW galaxy contains at least two sources, a main 
core source, and a fainter companion to the west/northwest 
at a velocity different by 600 km s^^ from the main NW 
core. This second galaxy in the interaction accounts for the 
unusually broad total CO linewidth (FWZP > 1500 km/s) 
of the NW galaxy merger. 

The far-IR luminosities of both sources are high, of 
order 10^"^ L0, indicating tremendous star- forming activ- 
ity, with estimated star formation rates greater than 1000 
M0yr~^. Since BR1202-0725 is radio quiet, its activity 
must be dominated by the massive starbursts that occur 
simultaneously in the galaxies interacting in the NW and 
SE components. We are obviously witnessing an extreme 
merging event in a group of galaxies, probably the youngest 
complex merger so far known. Future higher angular reso- 
lution observations of BR1202-0725 will allow us to better 
disentangle the properties of the individual starbursts that 
are interacting and merging, resulting in the remarkable 
appearance of this high-redshift object. 
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Appendix A: Absence of gravitational lensing 

Whether the submm/radio companion is a distinct object 
or a gravitational image of the quasar host galaxy is an im- 
portant question because in the absence of magnification, 
BR1202-0725 would be the most luminous binary CO and 
far-IR source in the Universe. The observational facts do 
indeed speak against gravitational lensing. The difference 
in hnewidths of the NW and SE sources, the absence of op- 
tical or UV emission from the NW source, and the muhiple 
CO and dust components in each source, are ah difhcuh to 
reconcile with a gravitational lensing scenario to account 
for the appearance of BR1202-0725. In the following dis- 
cussion of gravitational lensing, we definitely rule out the 
hypothesis that the SE and NW sources are multiple lensed 
images of the same background source. 

A.l. Different Hnewidths. 

The main piece of evidence comes from the greatly differ- 
ent molecular and atomic emission line profiles observed 
towards the NW and SE sources. These NW and SE CO 
and carbon-line profiles are so different in their widths and 
peak velocities, that they cannot be gravitational images 
of a single object. Differences in velocities and hnewidths 
might occur in multiple lensed images of optical lines, if, 
for example, a multiply-imaged quasar were locally mag- 
nified by a micro-lensing event that would affect only one 
of the images. Such an effect cannot affect the submm CO 
and carbon lines, however, because these low brightness, 
high luminosity lines are emitted over a large (> 2-kpc), 
relatively cool (50 to 100 K) region, and unlike the small, 
hot, accretion disk of a quasar, cannot be affected by micro- 
lensing. Differential image magnification between the NW 
and SE CO sources can also be ruled out; although the 
SE and NW CO sources have quite different hnewidths, 
their dust continuum fluxes are nearly the same. The dust 
sources and the molecular line sources are observed at the 
same wavelengths, and their measured sizes (~ 0.5") are 
roughly the same, so there is no evidence for any differ- 
ential magnification between the SE and NW sources at 
submm wavelengths. 
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A. 2. No optical counterpart of the NW galaxy. 

The second piece of evidence comes from the dissimilarity 
of the optical and mm/submm images of the quasar and its 
main companion, the NW source. The SE submm dust and 
molecule component is associated with an optically-bright 
quasar whose light is not seen at the position of the NW 
submm component. The nearly unsurmountable problem is 
to make a gravitational Icnsing model that would produce 
a single-spot optical image and a double-spot submm im- 
age. To first order, gravitational Icnsing is achromatic if the 
source size is the same at different wavelengths. 

But the quasar optical light (rest-frame UV) and its 
broad lines come from the black hole accretion disk on a 
scale of 10^^ pc, while the low-brightness-temperaturc CO 
lines and submm dust emission must be of kpc-scale size to 
be detectable at all at high redshifts. The molecular disk 
area is thus 10^° times larger than the quasar accretion disk 
area, and this can lead to dramatic differences in lensing 
amplification and lensed image shapes between optical and 
submm wavelengths. 

Two examples are IRAS F10214-F4724, and APM 
08279-1-5255. In F10214, the quasar is too obscured to be 
visible, but the AGN optical narrow-line region is magnified 
by a factor of 100, while the much larger submm molecu- 
lar disk is magnified by only a factor of 10. (see models by 
Broadhurst & Lehar 1995; Trentham 1995; Eisenhardt et al. 
1996). In APM 08279, the quasar's rest-frame UV through 
mid-IR is also magnified by factor of 100, while the cold part 
of the submm molecular disk is magnified by much less (see 
model by Egami et al. 2000, their Fig. 9). But even in these 
models with large differences in magnification between the 
optical and submm, it is nearly impossible to simultane- 
ously produce a double image at submm wavelengths, and 
only a single image at optical wavelengths. 

In a lensing scenario, one would have to imagine a 
galaxy on the line of sight somehow producing a double 
image only of the molecular disk, but not of the optical 
quasar, and then, in addition to this first lens that images 
only the molecular disk, a second lens prodiicing a strong 
micro- lensing of the quasar to make an optically bright, 
single, quasar image, with no micro-lcnsing amplification 
of the much larger molecular disk. This overly-complicated 
scenario would not work at near-IR wavelengths of 3.6 and 
4.6/im, (rest wavelengths 6300 and 8070 A), however, be- 
cause the emitting region, the hot outer parts of the cir- 
cumnuclear ring around the quasar, are large compared to 
a typical micro-lensing Einstein radius (e.g., Sluse et al. 
2011). 

Contradicting this scheme, the Spitzer archive images at 
3.6 and 4.6/U,m show that the quasar is indeed much brighter 
(>100 times) than any object toward the NW, in the near- 
IR. In fact, the light that is seen in the NW in Spitzer data 
is actually the extended starlight continuum from the Lya 1 
companion galaxy, which is also near-IR bright in the i?, /, 
and K bands (see discussion of their Spitzer data by Hines 
et al. 2006). The submm dust and molecular disks in the SE 
and NW galaxies that are discussed in the present paper 
cannot be detected in the Spitzer near-IR images, because 
they are too cold; their flux is hopelessly far down the Wien 
side of the blackbody curve, and they do not radiate at all 
at rest-frame optical wavelengths (Fig. 8). 

All these difficulties to find a lensing model that would 
produce a single-spot optical quasar image and a double- 



spot submm image, such as a galaxy-mass lens for a double- 
image molecular disk only, simultaneously with stellar-mass 
micro-lensing for the quasar only, that doesn't really work 
for the near-IR images, lead us to conclude that a lensing 
explanation is a blind alley. 

A.3. Absence of a candidate deflector to make a double 
image. 

Deep optical images from the HST and VLT do not show 
any evidence of a massive foreground object capable of de- 
flecting the light of the quasar of BR1202 0725 into two 
widely-separated images. Nor do they show any evidence 
of any arc-like structure, neither on galaxy-lensing angular 
scales, nor on cluster-lensing angular scales. 

Because the SE and NW sources have similar fluxes 
in the submm continuum and in the CO, [Cl], and [Cll] 
emission lines, the most natural hypothesis for gravitational 
Icnsing would be a fairly circularly symmetric mass distri- 
bution located close to the mid-point of the 4" line connect- 
ing the two sources. Deep ESO-NTT imaging (Fontana et 
al. 2000) of the system along with BVrIK ddXa. (Giallongo et 
al. 1998) allow us to put tight upper limits on the deflector 
brightness. 

To produce an image splitting of A6' = 4", a deflector at 
rcdshift zl, that for simplicity we approximate as a Singular 
Isothermal Sphere, would need to have a velocity dispersion 
given by: 

/ A 

a = 186.25 km/s ■ \ , (A.l) 

where w = d^s/ds is the ratio of the; angular diameter 
distances between the deflector and the source and between 
the observer and the source. We get a sensible value of tr < 
350 km s""'^ only if zions ^ 1- Furthermore, from the Faber- 
Jackson relation (e.g., Faber et al. 1987; Djorgovski & Davis 
1987), one can estimate the rest-frame luminosity of such a 
deflector and, for a SED typical of a massive galaxy, one can 
predict the redshifted r-band magnitude. For any redshift 
^icns ^ 1, the deflector would have a magnitude < 24 and be 
bright enough to be detectable in the HST image. We thus 
conclude that, unless the deflector is very atypical relative 
to more local lensing galaxies (e.g., Auger et al. 2010), one 
would be able to detect it. 

The brightest near-IR object in the vicinity of BR1202- 
0725 is the faint Lya 1 companion galaxy located between 
the NW and SE sources, which has an r magnitude of 24.3. 
It cannot play the role of a deflector, as it has the same 
redshift as BR1202-0725 and is much too faint. The second, 
even fainter, Lya 2 companion galaxy lies 3" southwest of 
BR1202-0725, but it is again at the same redshift as the 
quasar (Hu et al. 1996; 1997). 

Finally, there is no arc-like feature on deep optical im- 
ages that would indicate very strong magnification of the 
bright quasar. Such an arc-like feature might arise if the 
quasar (only) were on the infinite-magniflcation caustic of 
a lens so close to the line of sight to the quasar itself, 
that the Icnsing galaxy would be visible neither against 
the quasar light, nor even on images where the quasar light 
is subtracted out, such as the quasar-subtracted images of 
Fontana et al. (1996) and Hu et al. (1996). 

In summary, there is no evidence, at any wavelength, for 
a galaxy-mass lens along the line of sight to BR1202-0725, 
that could produce a submm double image. 
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A. 4. Absence of a nearby cluster lens for overall magnification 

Apart from trying to explain the submm double source by 
lensing by a galaxy-mass deflector, one might ask whether 
the entire BR1202-0725 complex could be lensed by fore- 
ground cluster of galaxies. In this case, the NW and SE 
submm sources are indeed independent galaxies, but each 
of them, and the optical quasar, might be amplified by 
a comparable amount, by a large, low-redshift foreground 
lens, i.e., a cluster of galaxies. Here, the answer is no. The 
BR1202-0725 region has been the object of a deep-field sur- 
vey, the NTT BR1202 field, together with an overlapping 
region, the NTT Deep Field, resulting in a total field of 
2.3' X 4' (Fontana et al. 2000). These two overlapping deep 
fields were surveyed in UBVr bands and in the near-IR 
IJHK bands. Magnitude limits were of order 26 in r and 
21.7 in K. Although there is a peak in the redshift distribu- 
tion at 2 ~ 0.6 (Giallongo et al. 1998), there is no evidence 
for a cluster of galax;ies that could act as an overall magni- 
fier for the BR1202-0725 complex as a whole. Nor is there 
any systematic lengthening of galaxy images into ellipses 
or arcs that would be signatures of large-scale magnifca- 
tion across the field. 
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